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Abstract 

We consider the generation of entropy when particle pairs are created at 
a cosmological level. Making a reduction via the particle number basis, we 
compute the classical limit for the entropy generation due to the evolution 
of the matter field fluctuations (squeeze transformation), obtaining that it is 
linear in the squeeze parameter for a general class of initial states. 

We also discuss the dependence of the generated entropy on the coarse 
graining criteria. 

PACS number: 98.80 

1 Introduction 

The cosmological particle production that occurs in early universe models, because 
of the changing space-time metric, could be a relevant source for entropy generation. 
In this respect, work has been done introducing an entropy measure that increases 
as the quantum field that represents these particles evolves. 

If we approximate the metric by a given time-dependent background, the evo- 
lution of the matter field fluctuations can be modelled; while, as noted by L. P. 
Grishchuk and Y. V. Sidorov , the corresponding Bogoliubov transformation that 
gives the relation between the in and out creation and annhilation operators can be 
expressed as a squeeze transformation. 

The generation of entropy due to particle production can also be studied in the 
formalism of squeezed states. The problem is that of obtaining in an unambiguous 
way, a notion of entropy sensitive to the unitary squeeze transformation S that the 
field undergoes. 

In Ref. H, B. L. Hu and D. Pavon associated the entropy generation with the 
increase of the mean particle number and the loss of coherence of the initial state. 

Alternatively, in Ref. ||, T. Prokopec considers the entropy as that obtained 
by coarse graining the density matrix, that is, by reducing the density matrix with 
respect to a given basis. In that work the entropy change is considered, when the 
initial state is the vacuum and the reduction is done via the occupation number 
basis or the (over-complete) basis of coherent states. In both cases, the classical 
limit (large squeezing limit) for entropy generation coincides. 

Another possibility, recently suggested by M. Gasperini and M. Giovannini 
is to use the basis of eigenstates of x, the superfluctuant quadrature of the field. In 
the x-basis the calculation can be exactly done, the generated entropy is linear in 
the squeeze parameter for all values of squeezing (see Refs. || and |J). 

In general, we could think in coarse graining a given density matrix p with 
respect to an observable A, having eigenfunctions \a), and define the entropy as 

S = — '^2(a\p\a)hi(a\p\a) (1) 

a 

Then, an important question that arises is to what extent the generation of entropy 
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due to the squeezing process depends on the initial state and the coarse graining 
criteria. 

In order to study these issues we will take the evolution of the system in the 
Schrodinger picture. If one considers a pure state \iJj), Eq. ([]]) reduces to 

s=-Ei«>i 2 H«>r (2) 

a 

This is indeed the Shannon entropy for the probability distribution, in the state 
associated with the observable A. The intuitive meaning of this expression is that e 5 
is the number of basis vectors which are "appreciably involved" in the representation 
of |V>) ("its richness"), see Ref. 0. 

Therefore, given a coarse graining criteria, we are interested in the way the 
richness of the initial state of the field changes, due to the parametric squeezing 
(production of particle pairs). 

In Ref. |I| each mode k of the field is independently associated with a squeeze 
transformation. A better description is given in Refs. and where a two 
mode transformation that mixes the modes k and —k is considered, conserving the 
momentum during the production of pairs. 

For the sake of simplicity we will first consider the change of the Shannon entropy 
under a one mode squeeze transformation \ip) — > S\ip), S = e~ lG , where 

G = U(a) 2 -a 2 ) , r > (3) 

(g^, a are the usual creation and annhilation operators). 

In Ref. [|] we studied the interplay between the quantum fluctuations in the 
superfluctuant (resp. squeezed) quadrature and the corresponding loss (resp. gain) 
of information. 

In this article we will study the classical limit for entropy generation showing 
that it is the same for a general class of initial states (here the coarse graining is 
done via N, the particle number operator). Secondly, we will study the dependence 
of the classical limit on coarse graining by considering some particular examples. 

In section §2 we review the eigenstates of G and their properties which lead to a 
very simple description of the state evolution. In section §3 we compute the leading 
order behavior (in the squeezing parameter r) for A = N considering a general 
class of normalizable initial states (which include particle number eigenstates). In 
section §4 we discuss the dependence of the entropy generation on the coarse graining 
criteria, and finally, in section §5 we compute the two mode case. 

2 The Eigenstates of the Squeeze Operator and 
the Propagation Kernel 

In order to study the physical consequences of the squeezing process it is convenient 
to use the wave function representation and in particular a set of eigenfunctions of 
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the squeeze generator (|3|). By defining the two canonically conjugate quadrature 
operators: 

^=y|( a + at ) > V=^jfi.^- a ) ( 4 ) 

we get a simple form for the squeeze generator G, that corresponds to a dilation 
transformation: 

r 

G=-(xp + px) , [x,p]=i (5) 

The eigenfunctions of the dilation operator @ have been introduced by C. G. Bollini 
and J. J. Giambiagi (see Ref. ||). By choosing the following realization for x and p 

1 d 

x = x , p = - — 

i ax 

the eigenvalue equation Gip = firip reduces to 

which leads to the solutions 

i)% = — =a:+~ 5 and ipt = -t=x_ 1 (7) 



where 

A _ J a; A if x > , a_/0 if a; > 

x +~) ifx<0 x_-< \x\ x iix<0 



Note that the spectrum is continuous and extends from — oo to +oo. 

The functions given in (0) form a complete set and satisfy ^-function normal- 
ization: 

dfiiiftixMlM + V^(xi)V_(x2)) = 6( Xl - x 2 ) (8) 

mm) =60* , (vciw = o (9) 

In Ref. we obtained the kernel associated with the squeezing operator S = e~ lG , 
i.e., K(x,x') = (x\S\x'), which can be rewritten using the completness relation: 

1 r+oo 1 .... . 

K{x, x') = ^~ d^^=(x^x' + - 1 ' 1 + x^x'S lli )e-^ r (10) 



oo \/X'X 



Therefore, if x' and x are greater than zero we obtain 

K(x,x') = — [ +0 ° dfi^=e i ^ nx - lnx '- r) 
27T J-oo yfx'x 



=5(lnx — lnx' — r) (11) 



' x'x 
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If x' and x are both negative, 



K(x,x') = -=<J(ln|x| -ln|x'| -r) (12) 



and if x' and x have opposite signs, we see from fllPp that K(x, x') is zero. These 
properties can be summarized by noting that the kernel K(x, x') is non zero only 
for x = x'e r : 

K(x, x') = e^5(x' - xe~ r ) (13) 

Analogously, we get: 

K(p,p') = e%5(p' -pe r ) (14) 

This can be directly seen from the symmetry of our system under the change x — > 
p, p — > — x, r — > — r (cf. (|)). 

From (p~3| ) we see that under squeezing, an initial state given by (x\tp) = ip(x) 
evolves into 

/+oo 
dx'K(x\x')i){x') 
-oo 

= e~2 / dx 5(x — xe~~ r )i})(x') 



Similarly, 



e 2ip(xe r ) (15) 
(p\S\ip) = e^if(pe r ) (16) 



where (p(p) = (p\ip). 

It is straightforward to verify that the normalization is preserved during the 
evolution. 

We can see that, for r > 0, x (resp. p) is the superfluctuant (resp. squeezed) 
quadrature: using (|l5|) and (|16|) we have for the dispersions in x and p, 

a r (x) = e r a (x) , a r (p) = e~ r a (p) (17) 

Note also that the eigenfunctions (0) are the analog to plane waves but instead of 
being invariant (up to a phase) under translations, they are invariant under dilations. 



3 Entropy Generation and the Initial State 

We will consider here the entropy generation which is obtained by coarse graining 
the density matrix. To do so we have to choose an observable A with eigenfunctions 
| a) and write the density matrix as 

p = EEI a >HpKX a 'l ( 18 ) 

a a 1 
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Then, this density matrix is reduced by setting to zero the off-diagonal terms to get: 

pred = ^2\a)(a\p\a)(a\ (19) 

a 

and the associated expression for the entropy is 

S = -Tr p red In p red 

= — y^(q|p|q) m(q|p|q) (20) 

a 

In contrast with the entropy computed with the full density matrix, which is unal- 
tered by squeezing, the entropy (^) is sensitive to this evolution and depends on the 
reduction scheme. A natural scheme is that associated with the occupation number 
basis, for which we will study the large squeezing regime (r >> 1) of the entropy 
generation. 

We will consider an initial pure state which is supposed to be normalizable. 
The corresponding density matrix p = leads to a Shannon entropy with 

respect to the particle number given by 

S = -£|»(r)>| 2 ln|<n|^(r)>| 2 (21) 

n 

where \tp(r)) = S\ip). Now, in order to obtain the large squeezing behavior of 
(PI), we first note that the amplitudes (n\ip(r)) go to zero in that limit, for in the 
^-representation we have (cf. (jT5|)): 

\(n\ip(r))\ = e 2 | / dxip n (x)ip(e r x)\ 



- 2 

hoo 



< e^M / rf J^llV (22) 

where M is an upper bound on the values of |-0(a;)|, ip n {x) are the eigenfunctions of 
the harmonic oscillator and H n (x) are the Hermite polynomials. As the integral in 
(P2|) is convergent, the overlapping (n\ip(r)) goes to zero at least exponentially for 
large r (and each term in (|2T| ) goes to zero at least as re _r ). 

Then, for large squeezing, the leading order of (21) comes from the infinite sum, 
and we can compute it by summing from a given no: 

oo 

S ~ - ]T |(n|^(r))| 2 ln|(n|^(r))| 2 , r»l (23) 

71=71(3 

Now, in order to make the computation of the asymptotic behavior easier we will 
consider the states in the p-representation, the quadrature which is being squeezed, 
and we will suppose that <p(p) is an even function (a similar reasoning applies for 
ip(p) odd). From Eq. ( |TED we have 



(2k\i/j(r)} = I dpip2k(p)ei(p(e r p) 

rS r+oo\ 



2e§ (/o + Is J d P^ipMe r p) (24) 
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where 6 is a fixed number. In the first integral we can use (see ||) 

f2k(p) = -^-r cos(2p^/fc +y 4 )(l + 0{±=)) (25) 

which is valid for large values of k (and p bounded). Therefore, 

(2A#(r)} = ~~~r (l +C ° ~ l +C °) dp cos(2p V ^+ I A)(l + 0(^=)Me r p) 
+2e* / dpv 2k {p)<p(e r p) (26) 



In order to set an upper bound on the integrals over the interval (5, oo) we demand 
that the associated wave function <p(p) satisfy, for large p, \<p(p)\ < const. p~ 2 and 
l^j^l < const. p~ 3 . In particular, this condition is fulfilled by the eigenstates of the 
particle number operator. We end up with an expression of the form 

<2£#(r)> = J|x J" dp co 8 (2pVk) (l + 0(-^)) V {e r p) + O (27) 

where we have replaced k + 1/4 by k within the same order of approximation con- 
sidered. Finally the change of variable p —>■ e~ r p leads to 

<2A#(r)) = f{ke- 2r )e- r + g{k e - 2r )0{e- 2r ) (28) 

with 

m = 2 r dp ^H v(p) (29) 

Jo y^jrki 

and g(k) being a function which for large k vanishes faster than 1/y/k. 

Now, replacing the sum by the integral, we can compute the leading order 
behavior of (1231): 

dk\(2k\i)(r))\ 2 ln\(2k\*P(r))\ 2 (30) 

Jko 



Choosing k so as to make valid the assumption fl25|), and using fl28|) , we see that 

/■oo 

S ~ - dke- 2r f 2 {ke- 2r ){-2r + \nf 2 {ke- 2r )) 

J kn 



~ 2r /°° dkf 2 (k)- I dk f 2 (k) In f 2 (k) (31) 

Jkoe~ 2r Jkoe~ 2r 

as fco is independent of r we can integrate the last expression from zero (in the large 
squeezing regime). From (p9|) we see that 



/■oo /-+oo ]_ /-CO . . - . 

/ dk f (k) = / dp' dp— / — cos (2v kp) cos (2v A;;/) <p(p')ip(p) 

Jo J-oo IT Jo hj V ' V 7 



,+oc ( I poo dk 

kl 

dp<p(p)<p(p) = 1 (32) 



The second term in ( |3~1"D is an approximation to the series for the initial entropy. 
Although this approximation is not so good, the difference with the exact series 
is r-independent and does not affect the leading order behavior of the generated 
entropy. Then, using (0), we finally get: 

AS = S - S ~ 2r , r » 1 (33) 

It is straightforward to generalize this result for a general initial mixture of normal- 
izable states \b): 

P = T,Pbb'\b)(b'\ . Trp=l (34) 

6,6' 

Following essentially the same steps that led from ( pT|) to (|3~3"D , for large values of 
squeezing, we also get 

S = -Y,Pbb'Y,( 2k \Kr))(b'(r)\2k) In \J2pw(2k\b{r))(b'{r)\2k) 

b,b' k \b,b' 

~ 2r (35) 

4 Entropy Generation and the Coarse Graining 
Criteria 

In the previous section we obtained that the generated entropy, when coarse grain- 
ing with respect to N, goes like AS ~ 2r, in the large squeezing regime. In general, 
if one considers another coarse graining basis, a different asymptotic behavior could 
be expected. For instance, if we take the singular case of the basis of the dila- 
tion operator eigenstates (cf. (0)), the associated entropy does not change under 
a squeeze transformation. This is precisely because these functions are stationary 
under squeezing. Now, let us take as an example the set of Laguerre functions 

1 _N 

Xn{x) = -j=e *L n (\x\) (36) 

which are similar (in the sense of their localization properties) to the harmonic os- 
cillator case, and let us consider the entropy change under a squeeze transformation 
of the state 

= ~ ¥ (37) 



dxe- ax LJx) = ± , 3?(a)>0 (3* 

o a n+i 



Using (see |0[) 



we have 

tan h r ' 



s 



and for the associated entropy we get 

S = -Y,\(Xnmr)}\ 2 \n\( Xn mr)}\ 2 

n 

= In (^cosh 2 (^ + sinh 2 (^) lntanh 2 (^) 

~ r (40) 

then, the leading order behavior of the generated entropy gives AS ~ r. 

There is however a factor of two with respect to the reduction done via the 
iV-basis. This can be qualitatively understood by using the intuitive meaning of the 
Shannon entropy of a state with respect to a given basis 0: e 5 is the number of 
basis vectors which are "appreciably involved" in the representation of As the 
squeezing developes, the width of the wave function in the ^-representation grows 
as e r . In the case of the iV-basis, in order to represent \if){r)), the \n) which are 



significant, can be estimated as those having J (x 2 ) n less than the width of the wave 
function we are considering. Taking into account that for the harmonic oscillator 



(x 2 ) n ~ «2, the number n of basis states that participate in the representation 

of \ip( r )) verifies ~ a^e 1 ' (a is the width of ip(x)) which gives S ~ Inn ~ 2r 
(cf. fl53|)). On the the other hand, in the case of the Laguerre functions, we have 
\J (x 2 ) n ~ \/6n, so applying the previous argument we obtain S ~ Inn ~ r (cf. 



We will now consider a coarse graining via the superfluctuant quadrature x. 
This is motivated by the approximated eigenf unctions used in Eq. (KSp which are 
plane waves (Dirac <5's in the x representation). This alternative basis was recently 
suggested by M. Gasperini and M. Giovannini. In this case, the calculation can 
be exactly done and the generated entropy turns out to be linear for all values of 
squeezing (see Refs. [[J and 0). In Ref. 0, this property was shown in the case 
where the inital density matrix is diagonal in the occupation number basis. From 
the state transformation given by Eq. ([H|), regardless the initial density matrix, the 
reduction done via the x operator leads to a linear generation of entropy under a 
squeeze transformation. If we make the reduction of ( j34|) using A = x in equation 
(|19D, we obtain: 

S = -Y,P»> j dx(x\b(r)}(b\r)\x)\n(j2Pbb'(x\b(r)}(b'(r)\x}) (41) 

6,6' \6,6' / 

where we have made explicit the squeezing of the initial state. Now, using ([15]) we 
see that 

(x\b(r))(b'(r)\x) = e- r Vb {e- r x)$ v {e- r x) (42) 
where <pb{x) is the wave function (x\b), so the expression for the entropy is: 



S = r I dx(x\SpS^\x) 

- I dxe~ r ^p W (p b (e~ r x)<pv(e~ r x)\nJ2Pbbi<Pb(e~ r x)Lpv(e~ r x) (43) 

J 6,6' 6,6' 
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The second term is just the initial entropy as one can see by changing variables 
e~ r x —>■ x; while the first term is rTr SpS^ = rTr p = r so we get: 

AS = r (44) 

This coincides with what we would have expected from our qualitative argument: 
the width of the state grows like e r , so the "number" of J's involved grows as e r , 
and the entropy goes like r. 

Finally, a reduction via the basis of coherent states can be considered. Starting 
from the vacuum, the leading order behavior can be calculated, giving AS ~ r. 
We can understand this behavior by noting that the localization properties of the 
over-complete basis of coherent states are similar to those of the Dirac delta's: the 
width is the same for every coherent state and they are uniformly distributed on the 
real axis. 



5 The Two Mode Squeeze Operator 

In Ref. m each mode k of the field is independently associated with a squeeze 
transformation. A better description is given in Refs. || and [|J where it is con- 
sidered a two mode transformation that mixes the modes k and —k, conserving the 
momentum during the production of pairs. 

The evolution operator can be expressed as (see [ |T^] ) 



U = RS (45) 



where 



R = e iJ dt>H (t>) ? H = n k (t)(aW + al k a- k + l) (46) 

fc,fcx>0 

S = e ld , G= £ ir k (t) (e- 2 ^a k a. k - e 2 ^a\al k ) (47) 

By defining new operators b k and c k (k x > 0) according to: 

a k = ^={b k -ic k )e~ ivar ^ , a ^ k = ±=(b k + ic k )e- i ^ (48) 

which satisfy the canonical commutation relations, S can be expressed in terms of 
two single mode squeeze transformations of the form given in @: 

S= II exp(^(6i 2 -6D)ex P (^(4 2 -4)) (49) 

Now, we will consider the generation of entropy under the evolution ( ^5] ) (in each 
pair of modes k and —k). For simplicity, we will first consider as coarse graining 
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basis: \nb, n c ) (the eigenstates of Nb = b\b k and N c = c\c k ) and then the more 



natural one, \n k ,n^ k ) (the eigenstates of N k = a\a k and N. 



a -k a -kj 



In the first case we have that the entropy at time t is given by 



S 



5^ (n 6 ,n c |t/" fc (t)|^) In (n 6 ,n c |£/fc(t)|^) 



(50) 



n b ,n c 



where U k {t) is the part of fl45|) that contains the modes and — fc. From 
and (P5| ) we see that the corresponding free part of the hamiltonian takes the form 
^k(t)(y k bk + c\c k + 1) which gives an irrelevant phase factor in the amplitudes of 

©■ 

Now, using an approximation similar to (|28|) , but for two single mode squezze 
operators, we obtain 



(n b ,n c \S k (tM) « e- 2r f(k b e- 2r ,k c e- 2r ) 
f(h, K) = 



OO rOO 



, cos(2p b y/h) cos(2p cA /fcJ 
dpb dp c 1 1 <piPb, Pc 



— OO J — OO 



(51) 



where n b = 2k b , n c = 2k c , S k (t) is the factor appearing in and we have sup- 
possed, as we did in (p8|), that the initial state is even. From (51]), we finally get 
AS(t) ~ 4r ft (f). 

The second case is more interesting because it corresponds to the occupation 
number basis, which has a direct physical significance. Let us consider for example an 
initial state having no particles with momentum k and no particles with momentum 
—k (the total momentum is zero). Again, the free part of (|45|) plays no role and the 
entropy at time t is 



S 



I A ^ /\ 

^ (n fc ,n_fe|Sfc(t)|n ,no) In (n fe , n_ fe |Sfc(*)|n , n ) 



(52) 



Taking into account that the squeeze operator S k (t) creates and annihilates pairs of 



particles in the modes k and —k (cf. pTf)), we have that the surviving terms in (52) 
are those having n k = n_ k = n: 



S = \{n,n\S k (t)\n ,n ) In (n,n\S k (t)\n ,n ) 

n 

We can change \n,n) to the basis \n a ,n b ): 
\n, n) = 



(53) 



. +CD" 1 " \/(2n - 2?')! (2?)! 



n!2 r ' 



i=o 



[n-j)\j\ 



where the prime in the second member denotes that the ket is given in the basis 
\n b ,n c ). 
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Using Stirling's formula we can approximate the factorials in as 



2 n (n-j)\j\ yft[{n-j)j]* 

moreover, for n large we shall replace the sum by an integral. Note however that in 
spite of the approximation we still have the normalization condition: 

* = h B( k- 5> = 1 • B ^ = TOT < 56) 

o Tr[(n-j)j]2 it 2 2 L[a + p) 

Now, for large squeezing, we obtain from (|5l|) 

(n ,n \Sl(t)\n,n) « e~ 2r f d 3 — -rf((n - j)e~ 2r , je~ 2r ) 

p -r ^ne~ 2r \ 

% Jo dj [(ne- 2r - j)j]^ (nC ~ J ' J) 
= e- r /(ne- 2r ) (57) 

Approximating the sum in (|53D by the integral and using (|57|) as well as the re- 
sult (see Appendix) J^°dnf(n) 2 = 1 (similar to (ft2])) we obtain the leading order 
behavior for the generated entropy in the occupation number basis: 

AS(t) ~ 2r fc (t) (58) 



6 Discussion 

In this paper we studied the dependence of the entropy generation on the initial 
state and the coarse graining criteria, relating these issues to the Shannon entropy; 
i.e., the change of the richness of a state in a given basis as the squeezing developes. 
For this purpose we considered the eigenstates of the squeeze operator which enabled 
us to describe the evolution as a dilation transformation of the wave function. 

Making a reduction via the particle number basis, we computed the classical 
limit for the entropy generation under a one mode squeeze transformation, obtaining 
that for a general class of initial states it is linear in the squeeze parameter. This 
is also the case when the reduction is done via the superfluctuant quadrature of the 
field (here, the linear behavior is valid for any squeezing). 

A deep understanding of entropy generation due to the creation of particle pairs 
will come from the knowledge of the physical mechanism that singles out a coarse 
graining basis. We note however that for a wide class of basis, the associated entropy 
will increase as the system evolves. This generation of entropy can be understood 
if we recall that the Shannon entropy gives the richness of a state with respect to 
a given basis. In the superfluctuant quadrature representation the wave function of 
the initial state flattens (cf. ([15])); and for a basis sharing the localization properties 
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of the occupation number basis or that of Dirac's (^-functions we will have that at 
each stage more basis vectors are necessary to describe the state. 

Then, we studied the entropy generation (in each pair of modes k and —k) when 
the matter field undergoes a two mode squeeze evolution, that is, when particle pairs 
are created or annhilated from the initial state, conserving momentum. Here, we 
considered the reduction via two different basis: the occupation number basis for 
the modes bk and Ck that factorize the two mode squeeze transformation, and the 
occupation number basis for particles with momentum k and —k. In the first case 
the classical behavior of the generated entropy is S ~ 4r^(t) (twice the value for the 
one mode case) while in the second case, starting from a state with definite particle 
number and zero momentum, it is S ~ 2rk(t). This calculation generalizes the result 
given in Ref. || where the initial state is the vacuum. 

Then, in the cases we have considered, we can see that the initial state does 
not leave its trace, when we look at the classical behavior of the generated entropy. 
This behavior turns out to be linear in the squeeze parameter with a proportionality 
factor that depends on the coarse graining criteria. 
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Appendix 

Here we shall show the result used in (|58|): 

dnf{nf = l (59) 



o 

In order to obtain an expression for f(n) we should find the equivalent of equa- 
tions (P5|) and (|29|) for the two mode case. 
First, using fl5"5|) and fl2"5|) we have 

1 /" ra 1 , 

(p b Pc\nn) « -r= / dj- —r(p b p c \2n-2j,2j) 

y/nJo [(n-j)j]? 

1 f n 1 / /r 

~ — / $ 77 TTT cos(2p b Jn - j) cos(2p c J j) 

ti-2 Jo [(n - 2 

^J (2^+^)) (60) 



Jo is a Bessel function. 

Now we can compute f(n) 



f(n) = J dp b dp c (n n \p b p c )-^=J (2yjn(pl + p 2 )j (61) 



1 no J(2n - 2j)U2j)\ r t , 

£ ^ / dp5dp c (p fe |2ra - 2jY(p c \2jyj [2^n{pl + pl 

~Tn \/7r n n — 7 !?! J V 



2"°^ AK-j)!j1 
1 ^ 1 /■ , , „ 



v% 1 /" Pg+Pg / / 

2 2n E tt^q- 7 ^ ^ C H2n »~ 2j ^ H2j (Pc)e 2 J °i 2 V H (Pb + PI 
/■oo 

= (-l) no / due- u/2 J (2y^)L no (u) (62) 
J 

where we changed into polar coordinates (p, a) in the plane p b ,p c and applied the 
identity (see fToH ) 

/ daH 2n (pcosa)H 2m (psma) = 2 7 r(-l)" + " A ! \ { 1 ' L n+m (p 2 ) (63) 
Jo n\m\ 

involving Hermite (H) and Laguerre (L) polynomials. We also used the normaliza- 
tion condition for the state \no,n Q ) in the basis \n b ,n c )' (cf. (|54f)): 

l_^(2n -2j)!(2j)! , 



With the above expression for f(n) we finally get: 

dnf{n) 2 = / dudu'e ~ / dn Jq(2^/tvu) Jo(2vra/)L no (-u)L no ('u') 
Jo Jo 



00 







due-Li, («) = 1 (64) 



where we have used the orthogonality of Bessel functions and the normalization of 
Laguerre polynomials. 
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